Magnetoresistance of atomic-sized contacts: an ab-initio study 
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The magnetoresistance (MR) effect in metallic atomic-sized contacts is studied theoretically by 
means of first-principle electronic structure calculations. We consider three-atom chains formed 
from Co, Cu, Si, and Al atoms suspended between semi-infinite Co leads. We employ the screened 
Korringa-Kohn-Rostoker Green's function method for the electronic structure calculation and eval- 
uate the conductance in the ballistic limit using the Landauer approach. The conductance through 
the constrictions reflects the spin-splitting of the Co bands and causes high MR ratios, up to 50%. 
The influence of the structural changes on the conductance is studied by considering different geo- 
metrical arrangements of atoms forming the chains. Our results show that the conductance through 
s-like states is robust against geometrical changes, whereas the transmission is strongly influenced 
by the atomic arrangement if p or d states contribute to the current. 

PACS numbers: 73.63.Rt, 73. 23. Ad, 75.47.Jn, 73.40.Cg 
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The investigation of electron transport through metal- 
lic atomic-sized contacts has attracted a lot of attention 
during the last 15 years, and most achievements are sum- 
marized in a recent review Q]. Using up-to-date ex- 
perimental techniques, such as mechanically controllable 
break junctions (MCBJ) or scanning tunneling mi- 
croscopy 0, 0], it is possible to fabricate nanocontacts 
with a single atom, atomic chain or a molecule 5] in the 
constriction. The experiments reveal step-like changes 
of the conductance upon elongation of the nanocon- 
tacts HUSH- In case of noble (Au, Ag, Cu) [SHU 
and alkali metals (Li, Na, K) 0, H, @ conductance his- 
tograms show a dominant peak very close to one conduc- 
tance quantum, G = 2e 2 /h, and smaller peaks close to 
integer conductance quanta. However, in case of transi- 
tion metals the situation differs significantly 0, 0] and 
a broad distribution of conductance values is usually ob- 
tained. 

To describe the nanocontacts theoretically, several 
methods have been developed by many groups during 
the last years. An approach based on the tight-binding 
(TB) Hamiltonian adapted for a nanocontact geometry 
was proposed in Refs. Il2t ll3L An important conclusion 
of TB models is that the conductance of single-atom con- 
tacts is related to the number of valence orbitals of the 
contact atom available at the Fermi energy [HE1. Lang 
and coworkers 0,0,0,0 and Kobayashi and cowork- 
ers 0, 0| studied the single- atom contac ts Fill , atomic 
chains of Al 0, Na [H 0, and C [Tjusing ab- 
initio calculations based on density functional (DF) the- 
ory with jcllium electrodes. The formation mechanisms 
of atomic chains made from different types of elements 
such as Ni, Pd, Pt, Cu, Ag and Au were investigated 
by means of molecular dynamics simulations [2lj . These 
studies were triggered by the experimental evidence of 
the formation of gold wires 141. Recently Mehrez et 
al. 0| and Brandbyge et al. |2flj presented fully self- 



consistent DF calculations of the conductance of atomic- 
sized contacts treating the electronic structure of both 
the contact and electrodes on the same footing. 

In this paper, we present ab-initio calculations of the 
conductance of nanocontacts and focus on magnetic sys- 
tems. We consider magnetic semi-infinite fee (001) Co 
electrodes joined by nanocontacts of different materials. 
We assume that they take the form of short three-atom 
Co-, Cu-, Si- or Al-chains suspended between the leads 
as shown in Fig. ^ Our aim is to investigate whether 
such hybrid systems could exhibit a large magnetoresis- 
tance (MR) effect. We wish to address the following the- 
oretical questions: (i) What is the influence of transition 
metal electrodes on the conductance of nanocontacts? 
(ii) What is the effect of the geometrical and electronic 
structure of the nanocontacts on the transport proper- 
ties? (iii) What are favorable conditions to increase the 
MR through a constriction? 

The structure of the nanocontacts studied below is 
presented in Fig. ^ The experimental lattice constant 
ao = 6.70 a.u. of fee Co was used in the calculations. 
In the first configuration (Fig. Hk), we consider a linear 
three-atom chain with the 1st and the 3rd atoms placed 
above the Co (001) surfaces at the ideal positions of the 
fee structure with the distance to the middle-atom of the 
chain being oo / y/2 which is the nearest neighbor distance 
in the fee lattice. In the second (zigzag-like) configura- 
tion (Fig. ^5)1 the atomic chain is a continuation of the 
fee structure along [001] direction, thus the symmetry 
is lower than in case of a linear chain. These two con- 
figurations were chosen to investigate the effect of the 
geometrical arrangement on the conductance. 

The electronic structure was calculated using the 
non-relativistic spin-polarized version of the screened 
Korringa-Kohn-Rostoker (KKR) Green's function 
method (for details, see Ref. E4I 1. The potentials were 
assumed to be spherically symmetric around each atom. 
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TABLE I: Magnetic moments (/Us) of atoms forming either linear or zigzag-like three-atom chains (see Fig. for parallel (P) 
and antiparallel (AP) orientation of the magnetic moments in the Co leads. The 'contact' atom of the chain sits above the 
(001) surface. The magnetic moment of bulk Co is 1.62 fiB, the moment of the Co (001) surface atom is 1.78 fiB- 



atoms 


Co linear Co zigzag 
P AP P AP 


Cu linear Cu zigzag 
P AP P AP 


Al linear Al zigzag 
P AP P AP 


Si linear Si zigzag 
P AP P AP 


contact 
central 


+1.91 +1.91 +2.00 +1.97 
+2.06 0.00 +2.28 0.00 


+0.08 +0.08 +0.08 +0.08 
+0.01 0.00 +0.03 0.00 


-0.09 -0.09 -0.10 -0.09 
0.00 0.00 0.00 0.00 


-0.11 -0.11 -0.07 -0.06 
+0.06 0.00 +0.11 0.00 



However, the multipole expansion of the charge density 
was taken into account. The angular momentum cut-off 
for the wavefunctions and the Green's function was 
chosen to be ^ max = 3 to ensure well converged results. 
Within the KKR methodic, the Green's function of 
the systems is obtained in two steps. First, we calculate 
the Green's function of the auxiliary system consisting 
of semi-infinite leads separated by a vacuum barrier. 
Next, the impurity problem is solved self-consistently 
by embedding the chain, surrounding atoms and empty 
sites into the host system, in order to account for 
the charge and spin density redistribution effects. We 
assumed parallel (P) and antiparallel (AP) magnetic 
configuration for the Co leads. The electronic structure 
of the constriction was calculated self-consistently for 
both cases. 

Ballistic conductance of the nanocontacts, g, was cal- 
culated using the Landauer theory as formulated by 
Baranger and Stone |2^|. We consider two semi- infinite 
Co (001) bulk leads connected by the scattering region. 
Conductance is calculated between two atomic planes lo- 
cated in the ideal leads. We neglect tunneling current far 
away from the contact region and sum up current contri- 
butions in real space in the vicinity of the constriction. 
Convergence was checked in order to achieve errors less 
than 5% in the conductance evaluation. A detailed de- 
scription and convergence properties of the method can 



(a) linear (b) zig-zag-like 

configuration configuration 




FIG. 1: Geometry of the nanocontacts: (a) linear configura- 
tion; (b) zigzag-like configuration. 



be found in Ref . H(| The real space version of the method 
employed here was presented also in Ref. 1271 

In Table I we present magnetic moments calculated for 
different atomic chains. We found that magnetic moment 
is enhanced for the Co atoms in the Co constriction reach- 
ing values up to 2.06 \xp for the central Co atom of the 
linear chain and 2.28 fip for the zigzag geometry. In the 
AP configuration, due to symmetry, the central atoms of 
the chains have zero moments and the magnetic profiles 
are antisymmetric. In case of Cu, Si, and Al the induced 
magnetic moments are small as can be seen from TablelH 

Our results for the MR at the Fermi energy (Ep) are 
summarized in Table IITI We define a magnetoresistance 
ratio as MR = (gp — gAp)/gAP where gp, gAP are the 
conductances for parallel and antiparallel magnetic con- 
figuration of the Co leads. The calculations predict a 
MR ratio of 30-40% for the Co constriction, and of 20% 
in case of the Cu chain suspended between the leads. The 
results for Al and Si chains are particularly interesting. 
We obtain MR around 50% for the linear configuration 
although magnetic moments are rather small (see Ta- 
ble H}. The values decrease to 20% when the symmetry 
of the constriction is reduced. 

In order to understand the relation between elec- 
tronic structure and transport properties we consider 
the energy-dependent transmission, T(E), through the 
constriction (Fig. [3 and Fig. 0). In the linear response 
the conductance per spin channel is related to the total 
transmission as g = e 2 /h J_™{—f'(E)}T(E)dE, where 
f'{E) is the derivative of the Fermi-Dirac distribution 
function. For zero temperature the conductance is g = 
e 2 /h T(Ep). However, in case of an applied voltage, V, 
states in the energy window of eV close Ep are relevant 
for the electron transport. 

In Fig.|21we present the transmission of Co, Cu, Al and 
Si constrictions between Co leads for both spin channels 



TABLE II: Magnetoresistance ratio MR — (gp~gAp)/gAP x 
100% for different zigzag and linear chains suspended between 
the Co leads. 
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FIG. 2: Spin-dependent transmission as a function of energy 
for different linear (left column) and zigzag-like (right col- 
umn) chains suspended between Co electrodes with magnetic 
moments oriented in parallel. The solid lines with up- and 
down-triangles correspond to transmission for spin-up and 
spin-down electrons, respectively. 



FIG. 3: Total transmission for both spins as a function of 
energy for different atomic chains suspended between the Co 
leads in case of parallel (black squares) and antiparallel (open 
squares) orientation of magnetic moments of Co. Transmis- 
sion for the parallel configuration is the sum of the transmis- 
sion for spin-up and spin-down electrons shown in Fig. [5] In 
case of the antiparallel configuration transmission is the same 
for both spins. 



for parallel magnetic configuration. We see that trans- 
mission for majority (spin-up) electrons is generally a 
rather smooth function of energy. However, for minor- 
ity (spin-down) electrons the transmission exhibits rather 
complicated behavior as a function of energy caused by 
Co d states which also contribute to electron transport. 
In Fig. [3] we present the total transmission (sum of both 
spins), for parallel (P) configuration in comparison with 
the antiparallel (AP) configuration. The conductance of 
the considered systems is the transmission at the Fermi 
energy, T(Ep), in units of e 2 /h. The difference of the 
transmission between P and AP configuration at Ep is a 
measure of magnetoresistance. 

Let us first concentrate on Co constrictions (Fig.[2^,b). 
The energy dependence of the transmission can be inter- 
preted in terms of local densities of states (LDOS). For 



this reason, in Fig. 0] we present the symmetry projected 
LDOS, s, p, d^ z 2_ r 2 1 d xz , d yz , at the Co central atom 
of the linear chain. In case of linear configuration the 
d xy , d x 2_ y 2 orbitals do not support current since they 
are oriented perpendicular to the current direction (z- 
axis) and form very sharp resonances in the LDOS of 
the central Co site of the chain because of weak coupling 
with the orbitals of the neighboring sites. From Fig. 0^ 
we see that Co majority states close to Ep are mainly 
sp-like since the d band is fully occupied and located be- 
low —0.75 eV with respect to the Fermi energy, while for 
minority states the Fermi level crosses the d band. This 
is also valid for the case of a Co zigzag chain. 

The examination of the LDOS at the Co sites of the lin- 
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FIG. 4: Symmetry projected local density of states (LDOS) 
at the Co atoms of the linear three-atom chain for the case 
of parallel orientation of magnetizations in the Co leads. We 
present the s, p and all d contributions relevant for a transport 
along the wire axis (z-axis). The contact Co atom is the chain 
atom next to the surface of the Co lead. Panel (a): s and p 
states of the central Co atom; panel (b): d 3r 2_ z 2 states of the 
central and contact Co atoms; panel (c): d xz and d yz states 
of the central and contact Co atoms. 



ear chain shows that due to the localized and directional 
character of the d orbitals in real space, transmission is 
high only if the orbitals of the same symmetry at neigh- 
boring sites of the chain are coupled. This can be seen 
from Fig. 0Ja and c. For example, the minority d xz and 
d yz states at — 1.2 eV (Fig.^J;) survive also at the neigh- 
boring atom (dotted line) and causes high transmission. 



Similarly, the jump in the transmission around —0.75 eV 
for spin-down electrons (Fig. |2Ji) correlate with the peak 
of the minority d 3z 2_ r 2 state (Fig. 0J)). 

Comparing Figs.|2t,b we can conclude that for Co con- 
strictions the transmission of spin-up electrons is rather 
insensitive to structural changes simulated by the linear 
and zigzag configurations of the atomic chains. On the 
contrary, the geometrical arrangement of the constriction 
seems to be more important for d electrons. Reducing 
the symmetry of the atomic chain by considering zigzag 
geometry destroys coupling between d orbitals, therefore 
the transmission of spin-down electrons is reduced and 
we obtain a smooth energy dependence (see Fig. Et>)- 

Our study is restricted to collinear magnetic configu- 
rations. To estimate the influence of this approximation 
we have considered also two-atom Co chains where a 180° 
domain wall is formed in the antiparallel configuration. 
In that case MR at E F was found to be approximately 
15% and the overall behavior is not far from the one dis- 
cussed above. Thus, our calculations predict a MR ratio 
rather sensitive to the geometrical and magnetic struc- 
ture. In addition, the spherical potential approximation 
might lead to small shifts of the electronic states to differ- 
ent energies, compared to a more accurate full-potential 
description. This might influence the MR values at the 
Fermi level reported in the paper. In particular, in cases 
where there is a strong variation of the transmission with 
energy close to the Fermi level, the reported MR values 
are expected to be less accurate. However, we believe 
that the configurations under consideration give a rep- 
resentative order of magnitude of the MR effect in Co 
nanocontacts, but cannot explain the MR ratio of Co 
contacts which was recently observed experimentally [13 • 
Moreover, our calculations for Ni nanocontacts using a 
zigzag configuration (Fig. ffls) give a magnetoresistance 
value of 24% contrary to the MR ratios of hundred per- 
cent observed in electrodeposited Ni nanocontacts. |2^| 
Thus, our calculations show that within the accuracy of 
density functional theory the experimental values cannot 
be explained assuming a clean, abrupt domain wall and 
ballistic transport. Further investigations in particular 
on the role of defects on magnetotransport are required, 
since a recent study indicates a significant enhancement 
of the MR ratio if oxygen is assumed to be at the contact 
region 29]. 

The next question is whether the electronic structure 
of the leads affects the transport properties through con- 
strictions? To answer this question the Co chains were 
replaced by Cu chains. It is well known that the majority 
band of Co matches the Cu band whereas the minority 
band of Co differs strongly from the Cu band. Our re- 
sults show that induced magnetic moments in the Cu 
chain are quite small (see Table [Q) and the spin-splitting 
of the LDOS is not significant (see Fig.EJ). Consequently, 
the calculated difference of the transmission for the two 
spins is mainly caused by the Co leads (see Fig. |2-\d). 
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FIG. 6: Symmetry projected local density of states (LDOS) 
at the Si atoms of the linear Si chain in case of parallel orien- 
tation of magnetizations in the Co leads. Solid line: central 
Si atom; dashed-dotted line: contact Si atom next to the Co 
surface. The peaks of the LDOS around —1.8 eV and near 
Ef have p x ,p y character. 



FIG. 5: Symmetry projected local density of states (LDOS) 
at the central Cu atom of the linear three-atom chain for the 
case of parallel orientation of magnetizations in the Co leads. 
Panel (a): s and p states; panel (b): d 3r 2_ z 2, d xz and d yz 
states. 



In case of the linear configuration there are highly con- 
ducting states around — 1.2 eV both for spin-up and spin- 
down channels. This is due to the matching of the Co 
d band to the Cu d states which are located in the en- 
ergy region below —0.75 eV (Fig. Eb). However, these 
highly conducting channels are closed in the zigzag con- 
figuration (Fig. I2J1) due to the same symmetry reasons 
discussed in case of the Co chains. Above —0.75 eV the 
electronic structure of spin-up Co states and Cu states is 
similar, thus the majority (spin-up) transmission of the 
Cu chains (Fig. [2c, d) is similar to the majority transmis- 
sion of the Co chains (Fig.|2Ja,b). The energy-dependent 
transmission for spin-up electrons is smooth in both con- 
figurations. However, the transmission of the spin-down 
electrons (Fig. |U,d) is significantly reduced and, finally, 
a MR ratio of ~ 15% at Ep is obtained (see Fig.^.d). 
The value 1.1 e 2 /h of the majority (spin-up) conductance 
at Ep for three-atom zigzag chain of Cu is in agreement 
with the conductance of an infinite monoatomic zigzag 
chain of fee Cu oriented along [001] direction. Since in 
that case only one band is crossing the Fermi level, the 
conductance is e 2 /h per spin channel. Finally, we con- 
clude that the electronic structure of the leads strongly 



influences the transmission through the constriction. The 
s-like states in the leads (majority states of Co) show 
a high transmission through an s-like chain, however, 
due to the band mismatch d-like minority Co states are 
strongly scattered and show a low transmission through 
the Cu chain. 

We further consider conductance through chains made 
from sp-elements connected to the leads. We have cho- 
sen Si and Al, since Si clusters of 2-20 atoms can be 
produced in the gas phase [3(|, thus measurements sim- 
ilar to Ref0 might be possible. Although the induced 
magnetic moments of Al and Si atoms are small, the MR 
ratio can be rather large (see Tables HllT)l . This is due 
to the " splitting" of the energy-dependent transmission 
curves which reflects the spin-splitting of the Co bands 
(see Figs. Fig.[2£-h). For example, in case of a linear 
Si chain (Fig. |2b) , the transmission of spin-up electrons 
increases up to 3 around —0.6 eV while the spin-down 
channel has low transmission which rises only above Ep. 
As a consequence, the spin-polarization of the current at 
Ep is high, P = (gj — gi)/{g] + g{) = 71%, which causes 
50% MR (Fig.Eb,h). The characteristic jump of - 2 in 
the transmission that is obtained in both spin channels at 
different energies can be explained by means of the spin- 
resolved LDOS at the Si sites (Fig. HJ). The majority 
LDOS of the central atom shows a pronounced peak at 
—0.6 eV below Ep shaped like a van-Hove singularity in 
one dimensional systems. This peak has p x ,Py character. 
The same p x , p y peak is seen in the minority band but it 
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is shifted due to induced spin-splitting to higher energies 
just above Ep. This state is responsible for the increase 
of the transmission by ~ 2. The minority LDOS of the 
central Si atom shows another pronounced resonance at 
— 1.8 eV below Ep which does not occur in the majority 
band and stems from the Co-Si interaction. This peak is 
also visible in the minority LDOS of the Si close to the 
Co surface. Since in the minority band the Si van-Hovc- 
like state becomes unoccupied due to spin-splitting, the 
electrons occupy the state formed by the Co-Si interac- 
tion which leads to nearly zero magnetic moments in the 
Si chain. The same effect of "splitting" of the energy- 
dependent transmission curves is also observed for Al, 
but the one-dimensional character of the p x ,p y state is 
less pronounced compared to Si, therefore the increase 
of the transmission of the linear Al chain is smoother 
(Fig. IB). 

In the zigzag configuration the p x ,Py degeneracy is 
lifted, so that for Si (Fig. |2h) the transmission increases 
in steps of ~ 1 both for spin- up and spin-down channels. 
In case of Al the steps are not so pronounced (Fig. Gf ) . 
comparing the zigzag with the linear chains we see that 
the " splitting" of the transmission curves is conserved in 
all cases (Fig.Et _ h). The MR ratio, however, drops for 
both Al and Si zigzag chains (Table 0. 

In conclusion, we have presented ab-initio calculations 
of the conductance of magnetic atomic-sized contacts. 
MR ratios up to 50% are predicted for different three- 
atom chains suspended between Co leads. Our results 
show that the conductance through s-like states is ro- 
bust against structural changes whereas the transmission 
is strongly influenced by the atomic arrangement in the 
chain if p or d states contribute to the current. Further- 
more, we have demonstrated that the spin-polarization 
of the ferromagnetic leads induces a spin-polarization of 
the current although the induced magnetic moments in 
the chain could be rather small (Cu, Al, Si). The in- 
duced spin-splitting of the wire states due to ferromag- 
netic leads is practically zero for s states but could be as 
large as 1 eV for p states. In case of sp chains (Al, Si) 
this causes strong spin-polarization of the current and, 
consequently, large MR values. 
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